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ABSTRACT 


In  the  portion  of  the  study  covered  in  this  report, 
linearized  theory  is  developed  for  potential  flow  around 
a  slender  body  at  not  too  large  angles  of  attack.  The 
solution  is  obtained  by  "definitizing"  the  flow  equation 
and  rigorously  satisfying  the  boundary  conditions.  This 
approach  is  completely  different  from  conventional 
methods  in  that  it  constitutes  a  process  of  obtaining 
an  analytic  solution  to  the  problem  of  potential  flow  by 
a  "march  from  the  body"  towards  infinity. 

The  extension  of  this  new  theory  to  the  nonlinear 
case  is  immediate.  Since  there  is  no  restriction  of 
the  Mach  number  made  in  Sections  5  and  6,  it  is 
expected  that  the  new  theory  should  be  applicable  to  the 
cases  of  subsonic,  transonic  and  supersonic  speeds 
as  well. 
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1.  INTRODUCTION 


The  three-dimensional  equations  of  motion  for  unsteady  isentropic 
flow  around  a  slender  body  are  linearized  under  the  assumption  of  angles 
of  attack. 

The  linearized  equations  of  four  independent  variables  are  reduced  to 
two  equations  under  the  assumption  of  a  potential.  These  two  equations 
relate  the  potential  to  the  condensation.  After  the  elimination  of  the  con¬ 
densation,  a  single  equation  for  the  potential  is  obtained.  This  equation  is 
called  the  indefinite  flow  equation.  It  is  valid  for  the  potential  around  any 
body  placed  in  the  stream. 

The  potential  is  written  as  a  sum  of  two  potentials,  the  steady  part  and 
the  unsteady  part.  The  unsteady  part  is  represented  as  a  product  of  a 
space  factor  and  an  oscillatory  time  factor.  The  general  approach  of  solv¬ 
ing  a  flow  problem  is  to  obtain  a  general  solution  of  the  indefinite  flow 
equation  first,  then  to  force  this  solution  to  fit  the  boundary  conditions  of  a 
body  under  consideration.  It  is  unfortunate  that  this  approach  has  not  been 
very  successful. 

Physically  speaking,  there  would  not  be  any  flow  problem  if  the  body 
were  not  placed  in  the  stream.  In  other  words,  the  flow  fields  are  created 
by  the  presence  of  the  body  through  its  definite  boundary  conditions.  The 
indefinite  flow  equation  is  too  general.  It  has  to  be  made  definite  for  a 
given  body.  This  process  is  called  the  principle  of  definitization. 

In  this  study,  the  principle  of  definitization  is  applied  at  the  outset. 

The  flow  equation  is  satisfied  near  the  surface  of  the  body  and  the  boundary 
conditions  are  taken  into  consideration  right  at  the  beginning.  The  solution 
thus  obtained  contains  an  unknown  factor,  which  when  introduced  into  the 
indefinite  flow  equation,  gives  rise  to  a  definitized  flow  equation. 

This  method  of  attacking  the  boundary  value  problem  for  the  potential 
can  easily  be  extended  to  other  boundary  v?.lue  problems.  It  is  completely 
new;  at  least  to  our  knowledge,  there  is  no  similar  approach  used  elsewhere 
in  the  literature. 

This  is  a  purely  theoretical  development.  No  calculation  is  made;  no 
application  is  attempted.  We  are  fully  aware  of  the  ways  in  which  the 
boundary  conditions  are  handled  and  manipulated  in  classical  potential  flow 
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theory,  and  of  the  facts  that  the  potential  solution  has  been  condemned  as 
"inadequate"  and  "unrealistic"  by  various  authors.  Consequently,  it  is 
hoped  that  this  new  approach  in  obtaining  a  new  potential  solution  might 
improve  the  classical  theory  and  restore  some  of  its  merits. 

This  report  contains  eight  sections.  Section  2  gives  the  fundamental 
equations  of  the  linearized  theory.  Section  3  derives  the  boundary  conditions 
for  unsteady  potential  flow  around  a  flexible  body.  Section  4  shows  the  way 
which  leads  to  arriving  at  the  new  approach  to  the  boundary  value  problem. 
Section  5  gives  the  solution  of  the  unsteady  potential,  Section  6  presents  a 
new  solution  to  the  classic  problem  of  steady  potential  and  Section  7  contains 
further  discussion  of  these  solutions. 

An  engineer  who  is  not  interested  in  theory,  may  read  Section  5  and 
Section  6  without  going  through  Sections  2,  3,  and  4;  however,  he  is  also 
cautioned  that  the  flexible  body  under  consideration  should  have  second 
order  partial  derivatives  with  respect  to  X£  and  X3,  and  that  the  assumption 
of  X!  =  g(x2»  X3)  as  the  body  equation  is  immaterial.  It  is  only  necessary 
that  one  of  the  .space  coordinates  can  be  expressed  as  a  function  of  the 
other  two.  In  case  of  corners,  special  treatment  has  to  be  made. 
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2.  FUNDAMENTAL  EQUATIONS 


PRESSURE  DENSITY  RELATION 

For  adiabatic,  isentropic  flow  the  pressure  density  relation  is 
given  by 

P  =  RpY  (1) 


where  R  and  y  are  constants. 

Let  pi  and  pj  be  the  free  stream  pressure  and  density  one  has  from 
Equation  1 


(2) 


One  may  set 


p  =  pj  (1  +  s)  (3) 

where  s  is  the  "condensation.  "  Substituting  Equation  3  into  Equation  2  it 
follows  if  s  is  small 

J-  =  (1  +  s}Y  =  1  +  Vs  (4) 

P1 


Consequently,  one  obtains 


1  Dp  P1V  D(1  +  s)  _  Dln(l  -f  s) 

p  Dt  -  p^(l  +  s)  Dt  p^  Dt 


Writing, 


(5) 


and 


In(l  +  s)  =  u4 


plV 


2 


a 

o 


(6) 


(7) 
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where  aQ  denotes  the  velocity  of  sound.  Equation  5  assumes  the  form 


1  Dp  2 

p  Dt  '  ao 


DU  4 
Dt 


NONLINEAR  EQUATIONS  OF  MOTION 


(8) 


Denoting  the  Cartesian  coordinates  by  x^ ,  x^,  x^,  the  time  by  t  and 
velocity  components  in  xj ,  X2,  x^  direction  by  Uj ,  U2,  U3  respectively,  one 
has  as  the  governing  equations  of  the  fluid  motion:  the  equation  of  continuity 


I 


k  =  1 


3  au\ 

2  Uk  9xk  ) 
k  -  1 


and  the  momentum  equations 


I 


k  =  1 


9  u.  _  9u. 

_ l  2  4 

Uk  9x,  +  ao  9x. 


0 


where  i  =  1,2,3 
LINEARIZATION 


(9) 


(10) 


Assuming  that  the  fluid  flows  over  a  body  (Figure  1)  at  an  angle  of 
attack  a  when  referred  to  a  certain  body  axis,  one  can  write 


u^  -  U  cos  a  +  u^ 

u2  *  a2  <U> 
u3  =  U  sin  a  +  u^ 

where  U  denotes  the  free  stream  velocity. 
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Substituting  Equation  11  into  Equation  9  and  assuming  that 


9u4 

9*k 


«< 


3u4 

"at" 


u  cos  a 


9xi 


(A) 


u  smft 


9u 

9x0 


4  f 


one  finds  by  dropping  higher -order  small  terms  the  linearized  continuity 
equation 


Figure  1.  Flow  Direction  and  Coordinate  System 
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Substituting  Equation  11  into  Equation  10  and  assuming  that  for 
i  =  1,  2,  3 


U  cos  a 


V  ~ 

1  uk  a xk 


k  =  1 


U  sin  a 


3u. 

2  _ 1 

Lo  ax. 


one  obtains  after  dropping  the  higher -order  small  terms  the  linearized 
momentum  equations 

9u.  9u.  au.  9u, 

i  „  x  .  i  2  4  „ 

— — -  +  U  cos  a  - —  +  U  sin  a  — —  +  a”*  — -  -  0 

at  ax,  8x„  o  ax. 

1  3  1 

i  =  1,  2,  3 

Equations  12  and  13  are  the  equations  to  be  investigated  in  this  report. 

EQUATION  FOR  THE  POTENTIAL 

Since  the  flow  is  irrotational,  one  can  assume  the  existence  of  a 
potential  <5  ,  such  that 


_  a  $ 

Uk  ”  3x 

k 


k  =  1,2,3 


Substituting  Equation  14  into  Equations  12  and  13,  one  obtains 


V  a2«  0U4  0U4  ®U4 

2,  — 2 - —  +Ucosa— -+Usi„a  — 

k  =  1  3xk  1  3 


1  /a $  „  a ®  .  9$ 

—  I— — +  U  cosff  ; - +  U  sin  a  — - 

c.  Vat  ax  axo/ 

i  1  3 


2 

a  u  . 
o  4 


*  i  =  1,  2,  3 
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Integrating  Equation  16  one  finds 


9  <£  „  3$  .  „  .  9$  2  2 

-r—  +  Ucosar  - - +  Usinar  — -  -  a  u  .  =  -  a  u  .  n  (17) 

9 1  9x^  9x^  o  4  o  4,0 

where  q  is  at  most  a  function  of  time. 

Solving  Equation  17  for  u^  and  substituting  the  result  into  Equation  15 
one  finds 


9Z$ 

3x? 


3Zg> 
3  x? 


Usinar 


9V  3> 
3x  33t 


+  Ucosar 


dZi 


8x18x3 


a  2  A 

+  Usinar - —  1  -  Ucosar 


3x 


l) 


2  2 

a  $  a  $ 

+  Jcosar - —  +  Usinar 


a2$ 


a  x^  at 


9x 


1 


dx  dx 
1  3> 


1(18) 


a2®  „  a2*  „  .  a2* 

- —  +  Ucosar  -x - —  +  Usinar 


at 


9x  ^ 9t 


2  3u410_ 

9x.3t/  a°  3t  ~  ° 


Equation  18  is  the  single  partial  differential  equation  for  the  total  potential, 
4>,  consisting  of  a  steady  part  and  an  unsteady  part. 

Since  one  is  mainly  interested  in  harmonic  oscillations,  it  is 
reasonable  to  set 

4>  =  i|/(xlfx2,x3)  +  <|>  (xlfx2,x3)  elu)t  +  a 2J*  u4  Qdt  +  At  +  B  (19) 
where  A  and  B  are  arbitrary  constants. 


W  riting 


U  U 

-  cosar  =  M  , - sins  =  M 

a  a  a  c 

o  o 


(20) 


one  finds  the  equation 


(1  .  mz)_JL^ 

a  „  2 


3  x 


1 


a  2 

2M  M  - - +  (1 

a  c  9  x^  9  x^ 


2  dZ\ii  32i 1/ 

M  ) - -  +  - — 


3  x. 


'  3x 


=  0  (21) 
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for  the  steady  part  of  the  potential  and  the  equation 


d-M  a2)J4- 


.  2  -  2 
w  3  d  ,,  w2,  3  d 

2M  M  — - - — +  (1  -  M  ) - — 

a  c  dxdx  *  c  a  2 

13  ox  ^ 


2 

+  (-f-)  d  =  0  •'  (22) 

o 

for  the  unsteady  part  of  the  potential. 

For  M  >  1  both  Equation  21  and  Equation  22  are  of  the  hyperbolic 

type. 
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3.  BOUNDARY  CONDITIONS 


SLENDER  BODY  AND  DEFORMATION 

Let  us  consider  a  flexible  slender  body  (Figure  2)  given  by  the 
equation 


£(x1'x2*X3)  =  °  (23) 

before  distortion,  where,  at  each  point,  P(x'},  X2>  X3)  on  the  surface  away 
from  a  neighbor  hood  near  the  nose,  x^  is  much  larger  compared  to  X£ 
and  x^. 


Figure  2.  Slender  Bcdy  Before  Distortion 

Now,  let  us  express  the  Carresian  coordinates  x^  and  X3  in  terms 
of  the  polar  coordinates  r  and  0.  Denoting  the  axial,  circumferential  and  the 
radial  displacements  Ax'p  rA0  and  Ar  by  u,  v,  w  (Figure  3)  and  the  coordi¬ 
nates  of  the  deflected  position  of  the  point  P  by  Xj,  X£,  x^,  one  has 


Figure  3,  Distortion  of  the  Cross  Section 

In  the  theory  of  shells,  it  is  found  that  u,  v,  w  are  functions  of  x^,  0, 
and  t.  For  harmonic  oscillation,  one  can  write 

—I  '  r>  \  iwt 

u  =  u(x  0)  e 

,  icot 
v  =  v(x  ,  0)  e 

w  =  w(x  1  ,0)  e^ut  (25) 

where  to  denotes  one  of  the  natural  frequencies  of  the  whole  system. 

Substituting  x^  ,  x^  ,  x^  into  the  original  equation  of  the  body,  one 

gets 

f(x  -  u,  x?  -  wcos0  +  vsin0,  x  -  wsin0  -  vcos0)  =  0  (26) 
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Using  Taylor*s  expansion,  it  follows  as  a  first  order  approximation  of 
Equation  26  the  following  result 


f(x1,x2,x3)  +  6(x1,x2#  x3)  elwt  =  0  (27) 

where  compared  to  f(x ^ ,  x2,  Xg)  the  magnitude  of  6(x  j »  x2 ’  X3^  *s  verY  small. 

The  function  6(x^,  x^,  Xg)  has  to  be  determined  either  from  the  theory 
of  shells,  if  such  a  theory  is  available,  or  from  tests. 

BOUNDARY  CONDITIONS 

Let  F(xj,  x2,  Xg)  be  a  function  of  three  variables,  x  ^x.^^x  and 
g(x  2 1  X3)  be  a  function  of  two  variables,  x2  and  Xg.  Let  us  assume  that 
F(g(  x2,x  3),  x2,  X3),  the  result  obtained  by  substituting  g(x2,  xg)  for  Xj  exists, 
then  we  will  denote  this  result  by  a  pair  of  brackets,  i.  e.  , 

[f|  =  F(g(x2!  x3),  x2,  x3)  (28) 


Taking  the  total  derivative  of  Equation  27,  one  gets  the  condition  of 
flow  tangency  to  the  deflected  surface  as 


for  xl,x2,x3  related  by  Equation  27. 

From  Equation  19,  Equation  29  becomes 


11 


^1  M  Line  - 

,x3j  A 


[J4_ 

0X3 


+  ico 


iwt 

e 


0 

(30) 


Multiplying  out  the  products  and  collecting  terms,  one  has 


Ucosa  - 


30j 

ax. 


3  f 
3x. 


8  ib 

a^T 


Usina  - 


8  ib 


3x 


icot 

+  e 


Ucosa 


86  | 

3  ib 

i 

36 

9x71 

lax  I 

+ 

lax0 

- 


Usina  - 


\ 


Hence,  we  have 


0 

(32) 


as  the  boundary  condition  for  the  steady  potential,  ib  , 


3f  | 

3  0 

3  f 

[  90  1  ,  [  9f 

3^  1 

3x  J 

3x  i 

*r 

3  X  2 

l 3 x  2  J  1 3x  3 

8V 

where 


F(x  2,  x3) 


(33) 


(34) 
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and 


36  | 

a  0 

1 

36 

30 

36 

ax,  1 

la  x. 

T 

lax  1 

lax. 

ax,l 

(35) 


as  the  boundary  conditions  for  the  unsteady  potential,  0.  Since  one  is  only 
concerned  with  potential  flow,  the  partial  derivatives  of  the  deflection  function 
6  have  to  be  very  small  so  that  the  condition  in  Equation  35  is  automatically 
satisfied  and  one  has  only  one  boundary  condition  for  the  unsteady  potential 
to  satisfy. 

EQUATION  AND  BOUNDARY  CONDITION  FOR  THE  STEADY  POTENTIAL 


To  summarize,  one  has  for  the  steady  potential,  d»,  flow  Equation  21 


(1  -  M2)  -2-4-  -  2M  M 
a  a  2 
a  ^ 


■  ~  _A_  +  (1  0  (2!) 

a  c  3  x„  3  x  .  c  „  2  _2 


3  1 


3x3  3x2 


and  the  boundary  condition  of  flow  tangency 


9f 


a  x. 


Ucosff  - 


a 

ax. 


af  i  [  a uj 
I  la  x. 


a  x. 


af 

dx . 


Usinff  - 


3  4* 

a  x„ 


-  0(32) 


where  the  brackets  are  evaluated  at  the  surface 


f(xi,x2*x3)  =  ° 


(36) 


which  is  a  first  order  approximation  of  f(xj,  x 2,  X3)  =  0 


Since  there  is  only  one  boundary  condition  at  the  surface  when  the 
condition 


Oj  =  0  (37) 

at  infinity  is  imposed,  Equation  21  will  have  a  unique  solution. 

EQUATION  AND  BOUNDARY  CONDITIONS  FOR  THE  UNSTEADY  POTENTIAL 
The  equation  for  the  space  factor  <f>  of  the  unsteady  potential  is 
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(22) 
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_9 _ <£_ 


3x 


2M  M  a 

a  c  3x„8x, 
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-2i  —  M 
a  c 
o 


+  (1 


d  0  *  3  0  /u\£ 

jt:  -2l— u*tt+  —  41  =0 

3  o  1  o 7 

The  boundary  conditions  for  cf>  are 

»  r  A  f  i  r  A  /A  1 

=  F(X2»X3) 


m2)±^+  Bf|_ 

8x3  8x  2 


f af  1 

f  df 

f  3  f  1 

f  I  . 

[  8f  1  f  80 

iaxj 

laxj 

4* 

|8x2I 

la*J  + 

I3X3]  t  3x3J 

(33) 


and 


3q  3  6 

3  0 

3  *  1  la  x  J 

13X, 

0 


(35) 


where  F  is  given  by  Equation  34. 

Mathematically,  Equation  22  with  two  boundary  conditions  (Equation  33) 
and  (Equation  35)  evaluated  at  the  surface  will  not  have  a  solution  if  some 
condition  is  imposed  after  the  body.  However,  one  can  consider  the  problem 
as  having  a  moving  boundary  so  that  Equation  33  should  be  satisfied  at  the 
fixed  surface  f(xj,  x2,  X3)  =  0  and  Equation  35  should  be  satisfied  at  the 
moving  surface  defined  by  Equation  27.  This  means  that  the  problem  has 
actually  only  a  single  boundary  condition  consisting  of  two  parts:  one  described 
by  Equation  33  and  the  other  given  by  Equation  35.  For  actual  computation, 
Equation  35  has  to  be  relaxed,  since  all  the  partial  derivatives  of  6  have  to  be 
small  to  allow  potential  flow  to  exist. 

TRAILING  EDGE  CONDITIONS 


It  will  be  shown  that  if  the  condition  in  Equation  35  is  relaxed,  the 
solution  cannot  be  unique.  To  uniquely  determine  the  solution,  one  will  have 
to  impose  the  condition  that  flow  be  tangent  to  the  surface  at  the  trailing  edge 
and  the  slope  of  the  streamline  is  continuous  in  the  region  of  the  trailing 
edge.  From  Equation  4,  if  the  subscripts -u,  and  l  are  used  to  denote  the 
upper  and  lower  surface  of  the  body  respectively,  one  will  have 


P1Y(SU. 


) 
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Kutta's  condition  says  that  after  the  body,  which  means  for  =  L  +  c  where 
e  is  an  arbitrary  small  number  and  g(X2»  x^)  ==  0,  one  should  have 


V  (Su. 


provided  that  s  is  small. 

From  Equation  17  and  6,  one  finds 


0 


8$ 

at 


+  Ucosa 


a$ 

a  x. 


+  Usina 


11L 

= 

j  ^ 

x. 


_  1  +  s 

?.  ,  u- 

a  l°g 


,  “  a  (s  -  s)  =  0 

1  +  s .  o  n>  l 


Substituting  $  given  by  Equation  19  into  the  left-hand  member  of  the 
above  equation  and  remembering  that  both  and  $  are  independent  of  t  one 
obtains 


3ijj  3^ 

cos  a  a„  +  sin  — 
oxj  9x ^ 

u 

+  eiwt 

t 

3  4*.  .  a  4>  .  w  , 

cc,soa5T;  +  Sln“ax3  +lu4' 

•  m 

CM 

O 

-<  |p 
+ 

u4»0  +  A 

**  C 

11 

o 

for  all  time  t.  Since  the  last  bracket  is  independent  of  position  it  vanishes 
identically.  Hence  one  concludes 


a  .  aiii 

cosor  — - +  smar 


3  x. 


3  x. 


a 


=  0 


(38) 


for  the  steady  potential  and 


30  .  30  .  w 

cosff  -  r' .  +  sincr  ■  +  l  —  0 

3x  3x3  U  y 


1L 


=  0 


(39) 


for  the  unsteady  potential. 
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4.  FORMAL  SOLUTION  TO  THE  BOUNDARY 
VALUE  PROBLEM  OF  THE  UNSTEADY 
POTENTIAL 


GENERAL  REMARK 


In  the  classical  theory  of  steady  potential  flow,  the  condition  of  flow 
tangency  at  the  surface  of  a  rigid  body  has  not  been  strictly  satisfied.  This 
is  true  whether  the  solution  is  obtained  by  using  the  method  of  source  or  sink 
or  of  doublets  or  by  applying  the  technique  of  characteristic  values. 

Furthermore,  the  constants  contained  in  all  classical  steady- state 
potential  solutions  are  actually  functions  of  location.  In  the  unsteady  case, 
it  is  expected  that  they  could  be  functions  of  time  as  well. 

To  obtain  a  solution  of  the  problem  of  unsteady  potential  satisfying  the 
boundary  condition  rigorously,  it  is  clear  that  the  conventional  method  cannot 
be  used  and  a  new  approach  has  to  be  tried.  In  this  Section,  we  will  indicate 
the  line  of  thought  which  leads  us  to  a  formal  solution  of  the  unsteady  potential 
problem,  while  the  explicit  solution  of  this  problem  will  be  given  later  in 
Section  5. 


REFORMULATION  OF  THE  BOUNDARY  CONDITIONS 


As  it  was  done  in  Section  3  (Boundary  Conditions),  we  will  denote  the 
result  of  letting  xj  =  g(x2»,  x3)  in  a  function  of  xj,  X£,  X3  by  using  the 
bracket  notation  [  ]  .  It  is  shown  that 


3  4 

=  4  x  (g,  x2,  x3), 

34 

3  x. 

3x„ 

L  3J 

'  3 

L  2 

x2, 


*  x  (g>  *2,  *3) 
1 


(40) 


Let  us  assume  for  a  moment  that  the  steady-state  potential  has  been 
uniquely  determined.  Solving  Equations  33  and  35  for  [3  4/3  x^]  and 
[  9  4/3xj1>  one  obtains  as  a  first  order  approximation 


3  4 

3  4 

9  x 

"  F3 

9x, 

3  J 

2 

3  4  ' 

=  F. 

3  <J) 

3  x 


+  L 


+  f. 


(41) 


(42) 
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where  F3,  £3  and  Fp  are  given  functions  of  X£  and  X3  alone.  Equations  41 
and  42  are  two  parts  of  the  boundary  condition  for  the  space  factor  of  the 
unsteady  potential  to  be  used  in  discussions  in  this  section. 

EXPRESSIONS  FOR  THE  SECOND  DERIVATIVES 

Using  the  bracket  notation,  one  has  the  following  identities  in  x^  an^  x3 
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(43f) 


Because  of  the  relation 
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Equations  43  a  through  43f  are  not  linearly  independent.  It  is  easy  to  show  that 
the  rank  of  the  system  consisting  of  Equations  43  through  43f  is  five.  It  can 
also  be  shown  that  one  can  eliminate  Equation  43a  and  consider  the  remaining 
five  equations  as  a  system  for  the  following  six  unknowns.. 
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<NJ 

CD 

1 _ 

r  ,2  1 

3  0 

2 

> 

3  x  9  x  _ 

» 

„  2 

3  x  .3  x  _ 
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3  x  3  x 

8xi  j 

1  £ 

L8X2  J 

2  3 

L3X3  J 

3  1 

To  obtain  a  sixth  equation  for  the  unknowns  in  Row  45,  one  solves 
Equation  23  for  x^  in  terms  of  X£  and  X3  and  makes  use  of  Equations  41  and 
42,  The  following  flow  equation  at  the  surface  of  the  body  is  obtained. 

•  VI 


(46) 


Equations  43  a,  b,  c,  d,  and  f  and  Equation  46  form  a  linearly  independent 
system  for  the  unknowns  listed  in  Row  45. 
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Solving  this  system  and  denoting 
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one  finds 
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8  20 
8  xJ 
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(53) 


where  the  A's,  B's,  C's  and  F's  are  given  functions  of  x^  and  x^. 

Equations  48  through  53  give  linear  relations  between  the  second  order 
brackets  and  the  unknown  function  C,  and  its  first  order  derivatives. 


COMPUTATION  OF  THE  HIGHER  ORDER  BRACKETS 


In  this  section,  how  to  obtain  the  third  order  brackets  is  demonstrated. 

The  fourth  and  higher  order  brackets  will  be  obtained  in  a  similar  manner. 
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•Introducing  the  unknown  function  £,  one  gets 
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where  the  A's,  B's,  C's,  E's  and  G's  are  given  functions  of  x ^  and  x^. 

Similar  linear  expressions  in  terms  of  E,  and  its  derivatives  up  to 
second  order  can  be  found  for  the  remaining  third  order  brackets  by 
differentiating 
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. 

with  respect  to  x^  and  x^.  There  are  12  third  order  brackets  connected  by 
12  linear  equations;  however,  these  equations  are  not  linearly  independent. 
Three  of  them  are  linear  combinations  of  the  other  nine.  To  find  a  non¬ 
singular  system,  one  augments  these  nine  equations  by  the  three  equations 
obtained  from  Equation  22  by  differentiation  with  respect  to  x^,  x ^  and  x 3, 
then  setting  xj  =  g  (x2,  X3).  This  makes  a  system  of  12  linearly  independ¬ 
ent  equations  for  t|ie  12  third  order  brackets.  This  system  is  sufficient  to 
determine  I  &  1 


3  x 


1 


The  fourth  and  higher  order  brackets  can  be  obtained  in  a  similar 
manner.  In  fact,  all  we  need  are  the  bracketed  derivatives 
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Assuming  that  all  these  brackets  are  computed,  then  one  is  led  to  the 
expression 


p(  Xj  ,  X^i  x^) 
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for  0  < 


<  L 
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where  L  is  the  length  of  the  body  as  a  formal  solution  to  the  problem  of 
unsteady  potential  flow  with  flow  Equation  22  and  boundary  condition  Equa¬ 
tions  33  and  35.  The  factor 

b  (xj,  x^»  x^)  5  Xj  •  g(x2>  x^)  (59) 

vanishes  identically  in  x?  and  x  only  when 

X1  =  g(x2’  x3} 


is  substituted  for  x^. 
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5.  DETERMINATION  OF  THE 
UNSTEADY  POTENTIAL 


The  argument  advanced  in  Section  4  has  led  us  to  a  formal  soLucion  of 
the  unsteady  potential.  To  actually  determine  this  potential,  one  part  of  the 
boundary  condition  has  to  Le  relaxed.  This  is  condition  Equation  35.  The 
best  justification  of  doing  this  is  that  the  partial  derivatives  of  6  have  to  be 
small  in  order  to  keep  the  flow  potential;  otherwise,  the  body  cannot  be 
smooth  and  the  flow  will  be  turbulent,  which  is  beyond  the  scope  of  our 
study. 

DETERMINATION  OF  THE  MEAN  SURFACE 

In  Sections  3  and  4  we  have  considered  a  body  given  by  Equation  23: 

£(X1‘  X2’  X3)  =  ° 

When  it  is  placed  in  an  unsteady  air  stream,  it  experiences  a  deflection 

6(xl>  x  ,  x  ) 

The  flexible  body  oscillating  harmonically  in  the  air  flow  is  described  by 
Equation  27: 


n  \  ri  \  ICOt 

£(xl*  x2-  X3)+  6(xl-  V  x3>  e  ‘  0 

Now,  let  us  suppose  that  we  are  interested  in  an  observation  period, 
Tq,  where  Tq  is  small.  One  procedure  of  obtaining  a  mean  surface  would 
be  to  take  the  time  average  of  Equation  27.  This  gives 


*2,  x  )  +  6(x  j , 


X3>I^ 


0  (60) 


Taking  the  real  part  of  the  left-hand  member  of  Equation  60,  one 

finds 


sin  u>  T 


£W(xr  X2'  X3>  s  f(xl’  X2'  x3)  + 


0 


wT 


n 


6(Xj,  x2,  x  )  0 
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as  the  mean  surface  of  the  oscillating  flexible  body.  Solving  the  above 
equation  for  x\  one  finds 


X1  =  ^X2:  X3^ 


(61) 


RECURSION  FORMULA  FOR  THE  COEFFICIENTS 

From  the  discussion  given  in  Section  4,  one  is  led  to  conclude  that  the 
space  factor  0  of  the  unsteady  potential  has  the  form 
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(62) 


for  0  ^  x^  <  L 

where  L  is  the  length  of  the  body 

The  coefficients  of  this  power  series  have  to  satisfy  the  following 
equations.  These  are  flow  Equation  22: 
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and  the  condition  of  flow  tangency  Equation  33: 
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where,  because  of  the  smallness  of  the  partial  derivatives  of  6,  F  is  set  to 
zero  and  condition  Equation  35  is  relaxed,  since  it  is  approximately 
satisfied.  The  brackets  indicate  the  fact  that  after  the  derivatives  are  taken, 
xj  is  replaced  by  g(x2,  X3). 
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Taking  the  derivatives  of  0,  defined  by  Equation  62,  one  obtains 
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On  substitution  of  Equations  62  and  63  into  Equation  22,  it  follows  by 
comparison  of  the  coefficients  of  (  -  g)n  the  recursion  formuLa 


(n  -*•  2)(n  +  1) 


1  -  M2  +  2M  M  + 

a  a  c  9  x  ^ 


2 


2(n  +  1) 


M  M  - -  + 

a  c  9  x^ 


(64) 


1 

2 


M 

c 


A 

n+ 


for  n  =  0,  1 ,  2,  .  .  . 


We  will  assume  that  the  function  g(x2»  x-j)  is  infinitely  differentiable. 
Equation  64  determines  An+2  when  An+i  and  An  are  known.  In  particular, 
it  defines  A2  in  terms  of  Aj  and  A0.  Consequently,  An+2  is  a  function 
of  Aj  and  AQ. 

DETERMINATION  OF  An 

Recursion  Equation  64  may  be  written  as 

(n  +  2)!  A  „  +  2b  (n  +  1)!  A  +  cn!  A  =0  (65) 

n+Z  n+1  n 

where 


b 


+  2M  M 
a  c 


(66) 


x 


M  M  - — 

a  c  9  x . 


9 


9  x 


2 


+ 


1 

2 


M 

a 
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and 


c 


2  3g 

M  +  2M  M  — + 
a  a  c  0 


(67) 


x 


are  Linear  functions  of  the  operators  3/3x2,  3/3x^,  9^/9x^  and 

3^/3x^  with  coefficients  which  are  functions  of  and  X3  •  Neither 

b  nor  c  contains  n. 

Setting 


A  -  u<n) 
n  r  (n  +  1 ) 


(68) 


one  gets,  on  substitution  of  Equation  68  into  Equation  65,  the  equation 


u(n  +  2)  +  2bu(n  +  1 )  +  cu(n)  =  C 


(69) 


So  far  as  n  is  concerned,  Equation  69  is  a  difference  equation  with  constant 
coefficients. 

To  integrate  Equation  69,  one  defines  an  operator  E  by 


E  u(n)  =  u(n  +  1 ) 


(70) 


Equation  69  can  be  written  as 


u(n  >-  2)  -  -  (2b  E  +  c)  u(n) 


(71) 


The  solution  of  Equation  71  is 

u(2m)  =  (-l)m  (2b  E  +  c)m  u(0)  (72) 


for  n  =  2m  and 


u( 2m  +  1 )  =  (-1  )r°  (2b  E  +c)mu(l)  (73) 

for  n  -  2m  +  1.  Substituting  Equations  72  and  73  into  Equation  68,  one 
obtains  expressions  for  An> 
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THE  UNSTEADY  POTENTIAL 

From  Equations  68,  72  and  73,  one  finds 
co 

<P  =  ^  (-Dm  r"1  (2m  +  DfXj  -  g)  (2b  E  +  c)m  u(0)  (74) 

m=0 

I 

co 

v-'  .  ,  .m  -1  2m  + 1  r  .m  . 

+  \  (-1)  T  (2m  +  2)(xi  -  g)  (2b  t  +  c)  u(l) 

m=0 


while  the  boundary  condition  given  by  Equation  33  yields 


’  af 

H  . 

3  f 

3  e 

a  f 

)A  -  1  8f 

a  a 

o 

3  f 

3  x 

3  x 

9  x. 

3  x 

3  x. 

)Ai  a  x 

3X  + 

3x 

2  J 

2 

L  3  J 

3 

L  lJ 

/  L  2J 

2 

L  3 

9  A 


a  x 


(75) 


or 


u(l)  =  hu(0) 


(76) 


where 


h  = 


3  f 


a  x 


2  J 


h. 

3  x„ 


a  f 


ax 


JLl 

3  x„ 


3  f 


3  x 


lj 


3  f 


ax 


2  J 


3  x. 


af 


a  x. 


ax 


(77) 


Hence ,  we  have 


co 


0  =  £  (-l)m  r-1  (2m  +  lKxj  -  g)C™  (2b  E  +  c)ri)  u(0)  (78) 

m=0 


,2m 


CO 


^  (-l)m  r  1  (2m  +  2){x[  -  g) 


2m+l 


(2b  E  +  c)m  hu(0) 


m=0 
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Symbolically  written,  we  have 


<p  =  cos 


(x1  -  g)(2b  E  +  c) 


1/2 


u(0> 


+  sin 


(xx  -  g)(2b E  +  c) 


1/2 


(2b  E 


TRAILING  EDGE  CONDITION 


Condition  39  for  the  unsteady  potential  can  be  written 


as 


cos  a 


3  <p 


1J 


+  sin  a 


3  <p 


=  0 


where  3  p  /3  xj  and  3  p  /3  x^  are  evaluated  at  =  L  +  e 

g( X£l  X3)  =  0  . 

From  Equation  62  one  finds 


3  <p 


3  x. 


00 


(n+  !)An+l  (L  +° 


n 


n=0 


3  <p 


3  x. 


00  3  A 

v- ■>  n  ,  T  .  n 

I  ^7" L  + ' 

n=0 


Consequently,  one  has 


3  A 


n 


(n  +  1)  cos  a  A  ,  +  sin  a  — — —  =  0 
n+1  3 x. 


or 


,  n 

A  =  (-l)n  T  (n  -*■  1)  (tan  or)n  -  A 

n  no 

ox. 


(79) 

c)  ^ ^  hu(0) 

approximately 

(39) 

and 


(80) 

(81) 
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particularly  for  n  =  1 


A 


1 


-  tan  a 


a 


A 

o 


(80a) 


and  for  n  =  2 


A 


2 


2 

tan  a 


A 

o 


(80b) 


From  Equations  76  and  77  one  finds 


9  f  * 

3Ao  / 

‘  8f  " 

-  tan  or 

'  9f  " 

9  x  0 

9x,  ( 

9  x  _ 

9  x. 

2j 

2  \ 

L  3j 

L  iJ 

/ 


9  A 
_ ( 

9  x. 


=  0 


the  solution  of  which  is 


A  =  C 
o 


H  ( x  ^ ,  x  ^ ) 


(81) 


(82) 


whers  H(X2i  X3)  is  a  known  function  uniquely  determined  by  the  shape  of 
the  body  and  C  is  an  arbitrary  function  of  H(x;?,  x}) 


Substituting  Equation  82  into  Equation.  65  for  n  =  0  one  gets  the 
following  equation  for  C 


a(x2’  xs)  C"  +  p(X2»  x3)C'+V(x2,  x3)C  =  0  (83) 

where  a,  p  and  y  are  functions  of  anci  x3  uniquely  determined  by 

the  geometry  of  the  body. 

Equation  83  can  be  solved  at  least  numerically  to  give 
C  |H(x2i  X3)|  .  It  is  clear  that,  in  general,  only  an  approximate  solution 
can  be  obtained  Further  discussion  on  this  solution  will  be  found  in 
Section  7.  Here,  we  point  out  the  fact  that,  as  a  result  of  a  physical 
problem,  Equation  83  should  have  at  least  a  characteristic  root  with  a 
negative  real  part  indicating  the  existence  of  a  region  of  stability. 
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6.  EVALUATION  OF  THE  STEADY  POTENTIAL 


EQUATIONS  FOR  THE  STEADY  POTENTIAL 


This  report  would  be  incomplete  if  the  steady  potential  is  not  reevalu¬ 
ated,  simply  for  the  following  two  reasons:  first,  the  classical  potential 
theory  does  not  rigorously  satisfy  the  boundary  condition  for  a  sharp-nosed 
slender  body  of  revolution  and  even  less  for  an  arbitrary  slender  body; 
second,  the  steady  potential  is  also  essential  for  the  unsteady  aerodynamics. 

For  the  steady  potential  •  one  has  flow  Equation  21 


(1  -  M  2)  .  2M  M  +  (1  -  M  2)  LjL  + 

a  2  a  c  ^  c  _  2  _  2 


=  0 


1 


3x39  x1 


x- 


3 


and  the  condition  of  flow  tangency  (Equation  32) 


where  the  brackets  are  evaluated  at  the  mean  surface  (Equation  61) 

X1  =  ®  ^x2'  x3^ 

RECURSION  FORMULA  FOR  THE  COEFFICIENTS 

Advancing  the  same  argument  given  in  Section  4,  one  can  also  assume 

CD 

4 >  =  ^  Bn  ( x2»  X3)  (Xj  -  g)n  for  0  <  Xj  <  L  (84) 

n  =  0 


where  L  is  the  length  of  the  body.  Differentiating  Equation  84  and  substituting 
the  results  into  Equation  21  one  obtains  the  recursion  formula 
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(n  +  2)  (n  +  1) 


(1  -  M  2)  +  2M  M  +  (1 

a  a  c  9 


M  2,  (3JLY  JlZ- 

c  l9xj  \?x2 


B 


n  +  2 


=  2(n  +  1) 


M  M  — -  +  (1  -  M  2)  +  7  f1  “  Mc2)“^_g2  *85) 

a  eg  x  '  c  9  x_  8  x  9x..  3x  2  c  a„  4 

3  3  3  3  £  o  x3 


1  9  g 

2  a  2 
8  x2 


B 


n  +  1 


f  2  2 

a  x3  a  x2 


B 


n 


Equations  85  are  sufficient  to  determine  B^,  B^,  .  . 

Using  the  boundary  condition  (Equation  32)  one  finds  also 


Bo- 


.  in  terms  of  B^  and 


where 


h*  = 


9f 


3x. 


af 

liL 

3x  _ 

3x0 

B  =  h*BQ  +  k# 


3f 


3x . 


ax. 


af 


3x. 


9x , 


is  an  operator  and 


k*  =  -U 


/|8£' 

9g 

\lax. 

3x, 

8  x  -  J  3  x 


3f 


a  x . 


af 


8x, 


cos  a  + 


(86) 


(87) 


I  sin  a)  (88) 

8x3I  J 


is  a  scalar. 

Through  Equations  85  and  86  all  the  B's  are  expressible  in  terms  of  BQ. 

DETERMINATION  OF  Bn 

Recursion  Equation  84  may  be  written  as 

(n  +  2)  !  B  .  +  2b*  (n  +  1)  !  B  +  c*n  !  B  =  0  (89) 

n  +  2  n  +  1  n 
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wh  -r e 


-1 


b*  =  - 


9g 


(1  -  Ma")  +  ZMUcf^-  +  (1  -  M  ) 
a  9  x„ 


2, (M . 
V8X3 


9g 

9  x,. 


(90) 


V.  M  — +■  ( 1  -  M  2)  — — 

a  c9x  U  c  ;  9x3  9x3  9x2  9x2 


and 


c*  = 


(1  -  Ma2)  +  2MaMc 


+  (1  -  Mc  ) 


±g_ 

,9x. 


i£ 

9x- 


x 


(91) 


(1  -  Mc2)  — 


9x. 


2  2  2  2 

are  linear  functions  of  the  operators  9 /9x£,  d/dx^,  9  /9x2  »  9  /9x^  with 
functions  of  x2  and  X3  as  coefficients.  It  is  assumed  that  the  function 
g(x2,  x^)  is  infinitely  differentiable.  Neither  b*  nor  c  *  in  Equations  90  and 
91  contains  n. 


Letting 


B  =  v(n)/r(n  +  1) 
n 

one  finds  from  Equation  89  the  following  equation  for  v(n). 

v(n  +  2)  +  2b*v(n  +  1)  +  c*v(n)  =  0 


Defining 


£  v(n)  =  v  (n  +  1 ) 

one  obtains  the  symbolic  solution  of  Equation  93  as 

v(2m)  =  (-l)m  (2b*  £+c*)mv(0) 


for  n  =  2m  and 

v(2m  +  1)  =  (-l)m  (Zb^  ^  +  c*)™  (^v(O)  +  k*) 

for  n  =  2m  -f  1 . 


(92) 

(93) 

(94) 

(95) 

(96) 
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Consequently,  one  has 


B2m  =  (-l)m  T~l  (2m  +  1)  (2b*  £  +  c*)m  BQ  (97) 

B,  (  =  (-l)m  T"1  (2m  +  2)  (2b*  E+  c*)m (h*B  +  k*) 

cm  +  1  (J 

as  expressions  for  the  B's. 

CONDITIONS  AT  INFINITY 

Substituting  Equation  97  into  Equation  84  one  has 

co 

Wxj,  x2,  x3)  =  J  (-l)m  r*1  (2m  +  l)(xj  -  g)2m  (2b*  £  +  c*)m  BQ  (98) 

m  =  0 

+  ^  (-l)m  r"1  (2m  +  2)  (xl  -  g)2m  +  1  (2b*  £  +  c*)m  (h*BQ  +  k*) 

m  =  0 


or  written  symbolically 

Mxj,  x-,,  x^)  =  cos  | -  g)  (2b*  £+  c*)l/2j  BQ  (99) 

+  sin  [(Xl  -  g)  (2b*  £  +  c*)1/2j(2b*  £  +  c*)" 1  /2(h*BQ  +  k*) 

Equation  99  represents  the  steady  potential  which  satisfies  the  flow 
Equation  21  and  the  condition  of  flow  tangency  Equation  32.  Besides  these 
equations,  the  steady  potential  i()  has  also  to  satisfy  the  trailing  edge 
condition  and  the  conditions  at  infinity.  The  conditions  at  infinity  are 

4*(  x  1 »  *2>  *3)  (i) 

bounded  for  all 


and 


L  <  x,  <  ® 


for 


^(xj,  x^>  x^)  0 


2  2 

x  c  .  x  L  —  co 
2  +  3 


(ii) 
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That  condition  (i)  is  satisfied  can  be  seen  from  Equation  99  by  substitut¬ 
ing  L  +  I  for  and  g  =  0.  To  satisfy  condition  (ii)  one  has  to  choose  in 
such  a  way  that 

(2b*  £  +  c*)m  Bq  -  0  (100) 

and 

(2b*  £  +  c*)“  (h*B  +  k*)  —  0  (101) 

2  2 

for  ail  positive  integers  m  as  x.  4-  x  —  co. 

Equation  101  is  only  possible  when  h*  has  bounded  derivatives  and 

(2b*  £  +  c*)m  k*  -0  (102) 

,  ,,  2  2 
for  all  positive  integers  m  as  x  +  x  —  <n  • 

To  >w  that  condition  in  Equation  102  can  be  satisfied,  let  us  take  the 
case  of  <  paraboloid  with  elliptic  cross  sections 

f  =  Xj  -  g(x 2 ■  x^)  =  Xj  - 

where  (3  and  V  are  fairly  large  constants. 

Here,  one  has 


_9g_ 
3x  . 


=  2 


_2  _9g_ 

p2'  8x3 


=  2 


-2  ’ 


9f 


0  X, 


=  1, 


0f 
9  x. 


,ig 

9  x. 


9f 
9  x. 


_9g_ 

3  x- 


Consequently,  we  have 


9  K 

U(cos  a  ~  sin  a) 

9x3 


9g 


9g 

,3  x. 


+  1 


x. 


=  U(cos  a  -  2 


sin  a) 


X,  x„ 

4  — +  4-— +  1 


It  is  obvious  in  this  case  that 


2  2 

k*  —  0  for  x2  +  x^  —  ® 
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As  a  second  example,  let  us  take  the  cone 


r  -  2  2  2  2 

f  =  x1  -  cot  6q(x2  +  X3  )  =  0 


where  0^  is  the  cone  hal£  angle.  Here,  we  have 


9f 


1  x . 


=  1, 


0f 


x. 


=  .  is- 


Sx. 


9f 


3x. 


ag 


Jg 

ax 


X  . 


-sr®-  =  cot  e 


0 


x2  +  X3 


_  Jg_ 
7  ’  9x3 


x . 


=  cot  6 


0 


/  2  2 
X2  +  *3 


/ 

V 


Consequently,  we  have 


x 


U(cos  a  -  cot  0 


0 


k*  = 


x2  +  x3 


sin  a) 


cot  0Q  +  1 


The  condition  k*  remains  small  can  only  be  satisfied  when  the  cone  half 
angle  is  small 

TRAILING  EDGE  CONDITION 

Condition  38  for  the  steady  potential  ijj  leads  to  the  result 


n 


B  =  (-l)n  r"1  (n  +  1)  (tan  af  — - B 

n  „  n  U 

6  x 


which  when  substituted  into  Equation  89  for  n  =  0,  gives  the  equation 


a  2b 


0 


a  2b 


q, - —  4.  a. 


0  *7  aBo  _  . 

+  r  a .  +  a-„B„  =  0 


T  a  2  2  9  x_a  x  ’  3  „  2  4  9  x  5  0 

a  x3  23  a  x  ^  3 


where  the  c's  are  given  functions  of  x2  and  x3> 


(103) 


(104) 
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Writing  Equation  103  for  n  =  1  and  substituting  the  result  into 
Equation  86  one  gets 


8  B 


0 


tan  a  ~ — —  +  h'’'B.  +  k'1'  =  0 
9x.  0 


(105) 


Combining  Equations  104  and  105,  one  has  a  system  of  two  equations  for  Bq 
which  can  only  be  solved  analytically  under  certain  conditions.  In  general, 
only  a  numerical  solution  can  be  obtained.  Further  discussion  on  this 
solution  is  continued  in  Section  7. 
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7.  UNIQUE  DETERMINATION  OF  THE  UNSTEADY  POTENTIAL 


In  Section  5  we  have  derived  a  linear  ordinary  differential  equation 
(Equation  83)  for  Aq  =  C,  the  solution  of  which  is  evidently  not  unique.  The 
theory  developed  thus  far  is  still  incomplete,  unless  one  can  uniquely  deter¬ 
mine  Aq.  This  will  be  done  in  this  section, 

ASSUMPTIONS  AND  DISCUSSION 

Let  us  consider  for  a  moment  x 3  and  X3  as  parameters  and  denote  the 
roots  of  the  characteristic  equation 


or  +  pr  +  y  =  0 

by  r^  and  r^.  Then  the  general  solution  of  Equation  83  is  given  by 

r^H  r^H 

C  =  A(x2>  x3)e  +B(x2,  x^)e 


(106) 


(107) 


Let  us  assume  that  for  xz^  +  X3^  —  co,  C— - 0  then  for  Aq  =  C  and  An  defined 
by  Equation  81  Equation  62  satisfies  all  the  equations  and  boundary  conditions 
imposed  on  it  so  far.  Since  A  and  B  are  arbitrary,  the  function  p  is  not 
unique. 


THE  CONDITION  OF  STAGNATION  PRESSURE 

To  uniquely  determine  the  unsteady  potential,  one  recalls  that  in  the 
classical  potential  theory  for  a  sharp-nosed  slender  body,  additional  condi¬ 
tions  at  the  nose  were  imposed  on  the  velocity  potential.  Since  these  condi¬ 
tions  are  not  available  for  a  blunt-nosed  body,  one  could  fix  the  arbitrary 
functions  A(  x2,  X3)  and  B(  x^,  X3}  by  predetermining  the  stagnation  pressure. 
Since  the  body  is  at  an  angle  of  attack,  the  coordinates  of  the  stagnation  point 
are  functions  of  a.  Let  us  denote  the  stagnation  values  by  the  subscripts  s 
and  assume  the  unsteady  stagnation  pressure  in  the  form 


s  i  wt 

7"  =  Pl+e  P2 
s 


(108) 
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From  the  nonlinear  Bernoulli's  equation  one  has 


iiot  1  2  3  vp  tt  3  ^ 

P,  +  p_.e  =  —  U  +  U  cos  a- -  +  U  sin  a  - — 

1  2  2  9  Xj  9  x, 


+  e 


iwt 


tt  3  9  tt  •  3  <p 

l  9  +  U  cos  a  - — —  +  U  sin  a  — 


3  x 


1 


3  x 


3 


(109) 


By  comparison  one  has 


and 


1  TT2  TT  0^ 

~  U  +  U  cos  a  — — 

2  3  x 


+  U  sin  a 


3  i|j 


1 


=  iuj9  4-  U  cos  a 


3  9 

3  x. 


+  U  sin  a 


3  9 

FT 


(110) 


(ill) 


Since  pj  and  p £  are  given,  Equations  110  and  111  are  sufficient  to  determine 
the  arbitrary  "constants"  contained  in  t|j  and  9.  This  finally  leads  to  a  unique 
solution  for  the  problem  of  unsteady  potential  flow  around  a  slender  body  at 
angles  of  attack. 
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8.  CONCLUSION  AND  RECOMMENDATION 


The  results  given  in  this  report  represent  an  attempt  to  solve  the  prob¬ 
lem  of  unsteady  potential  flow  around  a  slender  body  first  by  satisfying  the 
boundary  condition  at  the  body,  then  by  working  the  solution  toward  infinity. 
With  regard  to  the  problem  of  convergence  of  the  solution,  no  positive  state¬ 
ment  can  be  made  at  this  time. 

To  evaluate  the  merit  or  demerit  of  this  method  of  solution,  it  is 
recommended  that  it  be  tried,  for  example,  on  a  paraboloidal  body  with 
elliptic  cross  sections  to  see  whether  or  not  the  series  converges.  If  the 
series  does  converge,  then  this  method  should  be  explored  further,  and  a 
convergence  proof  should  be  established  by  imposing  appropriate  conditions 
on  the  function  g  (x2»  X3).  If  the  series  diverges  for  a  paraboloidal  body, 
then  this  method  should  be  disregarded. 
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